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Abstract. Using the theory of optimal rocket trajectories in general relativ- 
ity, recently developed in IHN11I , we show that the "obvious" manoeuvre of 
using a tangential instantaneous acceleration to escape a stable circular orbit 
in the Schwarzschild spacetime satisfies the optimality conditions if and only 
if the magnitude of the acceleration is smaller than a certain bound. 



Introduction 

Using Newtonian mechanics, the energy per unit mass acquired by a rocket on 
an instantaneous acceleration which changes its velocity from v to v + Av is 

A£ = \ \v + Av\ 2 - ±\v\ 2 = ±\Av\ 2 + (v, Av). 

Therefore, the energy is maximized when the instantaneous acceleration is parallel 
to the initial velocity. 

Now assume that the rocket is on a circular orbit of radius r around a spherically 
symmetric body of mass M. Its initial speed is theifl 

An obvious strategy to escape the gravitational field of the mass M is to accelerate 
instantaneously in the tangential direction, thus maximizing the energy per unit 
mass. However, as first noted by Lawden [Law53 , this can only be the optimal 
strategy if \Av | < V. (The range of possible values of \Av | is a characteristic of the 
rocket, dictated by its fuel supply and ejection speed). This is easy to understand 
in the case when the central body is a point mass: indeed, if the rocket is capable 
of \Av\ > V, then it can instead use part of its fuel to reduce its velocity to an 
arbitrarily small value, entering an elliptical orbit whose pericenter is arbitrarily 
close to M; its velocity at the pericenter is then arbitrarily large, and a second 
tangential acceleration results in an arbitrarily large energy gain. 

Noticing that the rocket's angular momentum (per unit mass) after the instan- 
taneous acceleration is L = r{V + \Av |), we can express Lawden's bound as 
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In this paper, we use the theory of optimal rocket trajectories in general relativity, 
recently developed in |HNllj . to obtain the general relativistic version of |T]): 

(2) — < Hi: Hi 

r* - (i_3M) 3 _|(i_6M)(i_2M) 2 • 

This bound assumes that the rocket is initially on a stable circular orbit (r > 6M). 
Note that ((2]) reduces to (p} in the nonrelativistic limit y C 1; the two formulas 
differ significantly for circular orbits close to a compact central object, like a neutron 
star or a black hole. Due to the more complicated structure of general relativity, 
there is no simple interpretation of ([2]), even in the case when the central body is 
a black hole (the relativistic analogue of a point mass). 
For comparison's sake, we must have 

r 2 2M 

(3) ^> 



2 — f-i 2M 



to escape the gravitational field. Since the right-hand sides of $2$ and © become 
equal for r = (5 + \/7)M, we see that for r < (5 + \/7)M the optimal manoeuvre to 
escape the gravitational field is never a single tangential instantaneous acceleration. 

The structure of the paper is as follows. In Section[T]we summarize the theory of 
optimal rocket trajectories in general relativity, for the convenience of the reader. 
In Section [5] we obtain the endpoint conditions for the problem of maximizing the 
rocket's energy in a stationary spacetime. Finally, the bound © is deduced in 
Section [31 



1. The rocket problem in general relativity 

In |HNllj . the classical Newtonian theory of optimal rocket trajectories, devel- 
oped in [Law63j . was generalized to the general relativity setting. There it was 
shown that optimal trajectories are continuous, sectionally smooth timelike curves, 
obtained by piecing together free-fall (geodesic) and accelerated arcs, possibly with 
instantaneous (Dirac delta) accelerations at the junction points. They must satisfy 
the differential equations 







P 




= ~t + 







( 4 ) \ VuPi 1 = -q» + paU^ 



where U is the four- velocity, a is the magnitude of the proper acceleration, P is an 
auxiliary vector field (called the primer), p is the magnitude of the primer, and 
Rafi^v are the components of the Riemann curvature tensor. All quantities above 
are continuous at a junction point except possibly U and P, which change by a 
boost with positive parameter in their common plane if there is an instantaneous 
acceleration. The magnitude of p the primer, however, is C 2 at the junction points, 
and satisfies p < po on free-fall arcs and p = po on accelerated arcs and instan- 
taneous accelerations (for some constant p a > 0). Moreover, the equations above 
admit the first integrals 

(5) P^ - - 0. 
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The primer P and the vector q are related to the momenta m M and p^ conjugated 
to the variables x^ and by 

(6) P^P^ + {p a U a )U^ 
and 

(7) q^^m^-T^paU 13 , 

where rjg are the Christoffel symbols. If the final values x± and t/f of the variables 
x M and are not specified, then the final values {m f j,) 1 and {p^) 1 of the momenta 
must satisfy 

dJ_ 

and 

( 9 ) Om)i = -^> 

where J = J(x°, . . . , a;f , U®, . . . , t/f) is the quantity to be minimized. We shall 
make use of these formulas in the next section. 



(8) (%)i = -; v, , 



2. Maximizing the final energy 

To be able to speak of energy we consider a stationary (d + l)-dimensional 
spacetime with future-pointing timelike Killing vector field T. The rocket's energy 
per unit mass is then the quantity E = —T^U^, and in order to maximize its final 
value we minimize 

J(xi, . . . , Xi , U\ , . . . , f/j ) = T^x-y, . . . , Xi)Ui . 
If we do not fix x± and C/f then ([8]) and ((9]) imply 

3 J 

and 

(Pm)i = = — (-^")i' 

or, in view of and (J7J, 

(10) (^)i = - (r a c/ Q )^) 1 

and 

(11) (^) 1 = (-(^r,)c/^ + r^r Q ^) 1 

= (-(v M T„)EHi = ((v,r M )c/-) 1 = (v^, 

where we used the Killing equation V^T^ + V^T), = 0. 

Since energy is conserved along geodesies, we can always assume that the optimal 
trajectory ends on a free-fall arc. Along such arcs we have a = 0, and so ((4]) reduce 
to the geodesic and Jacobi (geodesic deviation) equations: 

VuU^ = 
VuP^ = -q>> 
V v % = R, iaf)l U a pPW 
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We can obtain a solution to these equations by choosing 

\Py. = -T„ - {T a U a )U^ 
1 Qp. = V[/T M 



(12) 



as T a U a is constant along the geodesic and both T and U are trivially Jacobi fields. 
Since these choices also satisfy the endpoint conditions (ITU)) and (fTTj) , they must be 
the solution in the final free fall arc. 



3. Optimality condition 



We now focus on the problem of escaping a stable circular orbit in the Schwarzschild 
spacetime. For simplicity, we restrict ourselves to motions on the (2+l)-dimensional 
totally geodesic submanifold corresponding to the equatorial plane, whose metric 
in the standard Schwarzschild coordinates (t, r, ip) is 



ds 2 = - 1 - 



2M 



dt 2 



dr 2 + r 2 dip 2 



The geodesic equations in this coordinate system can be written as 

-l 



(13) 



t + 



2M 



1 - 



2M 



2M\ . 

//• = u <= | 1 j t = E 



r 2 \ r J r 2 \ r ) 

2 2 
Cp H — rip = <^ r ip = L 



f _ r 



2M\ . 2 
1 ) <p 



r J 



where E and L are integration constants (energy and angular momentum per unit 
mass). Using the normalization condition 



(14) 



2M 



2M 



- 1 * + 1 r + TV = -1 



we find for a circular orbit traversed in the positive sense (f = r = 0, (p > 0) 



(15) 



2M 



3 M 



L = 



' Mr 

I 3M 



The Jacobi equation for such a circular orbit was studied in |HNllj . There it was 
found that the radial and tangential components of the primer satisfy 



(16) 



P' 



r J r 1 



2L 
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and so, for stable circular orbits (r > 6M) we have 

' P r = Asin(W) 

or A 

= -^-cos{ujt) + B 

r 6 u 

where A and B are integration constants, r is the proper time (chosen such that 
P r = for t = 0), and 



6M 
r 



(17) uj = — \/]- 

It was also shown in [HNllj that the magnitude p of the primer satisfies 



2= ( 1 _2M\ 1 



Therefore we find for r = 
(18) p 2 - 

p = and 



r ) 

) 

2M 
r 



(pry 



1 - 



3A/ 



2LA 



{P v y 



B 



(19) 



pp - 



r J 



A 2 u 2 - l - 



3M\ 2LAuj (2LA 



— - + B 



From (fI5]l, ([TT1) and (0, and assuming P v > for r = 0, we find 



2M\~ 5 / 3MV 
r 



and so 

(20) 

where 



2A 



r I 



'1- 



3M 
r 



Br 



(P-D), 



D 



2M 

r 



Using (HHJ) and ([50^. we can write as 
(21) PP 



= (l-™) " AulAu-(l- —) Vl - ' 



r / 

(1 



(/>-£>) 



— 



Now consider the manoeuvre in which the rocket escapes the circular orbit 
through a tangent instantaneous acceleration in the positive sense. Since for such 
a manoeuvre to be optimal we must have P r — and P v > (recall that the 
primer determines the direction of the instantaneous acceleration), we can assume 
without loss of generality that the instantaneous acceleration takes place at r = 0. 
For t > we then have, by (|T2"j) . 



dt 



EU, 
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where E is the energy per unit mass of the final free-fall arc, and so 



(22) 



P 2 = PaP" 



r J 



E z 



2E 2 = E 2 - 1 



2M 
r 



Since f = for r = 0, we have 



(23) 



M . 



PP 



AP 



2M\ 1 , 



r 

where we used (| 13[) . Since by (|T4")) we have 

2M 



1 



E 



AI 



we can rewrite 

(24) 

and 

(25) 



and (j2"3")l as 



2 _ / 2A1\ L 2 

I r J r 2 



M 



r / r 



Since p must be C 2 , the values of pp from expressions (|2ip and ([23)1 must be the 
same. Equating them and using we obtain 



(26) 



(1 



3M \ (-1 _ GAI \ 
Li„2 , St r_±< _ D) 

(l — 2M.y h — M£y ' 



4(1 -hL (p~ d )-p 



M 
r 



Notice that p must have a local maximum at r = 0, and so we must have p < 0, 
which by (|21[) implies ^4 > 0, that is, p > D. On the other hand, the maximum 
must be the global maximum of p over the circular orbit, since if one waits for an 
arbitrary amount of proper time before performing the instantaneous acceleration 
the trajectory thus obtained must still be optimal. Therefore we must have B > 0, 
that is, D > 0. 

Consider the function 



We have 



f(D) = {p-D) 



f(D) 



1- 



6M N 



4(l-% (p_jD) - p 



6M» 



2(1 



3M 



■(P-D) 



and so /'(£>) > for p > 0, r > 6M and D > 0. Consequently /(£>) > /(0), and 
so (|26| implies 



(1_3M) (1 _6M) 



(1-gM) 
4(1-%) 



< 



M 



which, using (p4|) . is easily seen to be equivalent to 



OPTIMAL ESCAPE FROM CIRCULAR ORBITS AROUND BLACK HOLES 



7 



References 

[HNll] P. Henriques and J. Natario, The rocket problem in general relativity, J. Optim. Theory 

Appl. 154 (2011), 500-52. 
[Law53] D. Lawden, Escape to infinity from circular orbits, J. Brit. Interplan. Soc. 12 (1953), 

68-71. 

[Law63] , Optimal trajectories for space navigation, Butterworths, 1963. 

[MTW73] C. Misner, K. Thome, and J. A. Wheeler, Gravitation, Freeman, 1973. 

Centro de Analise Matematica, Geometria e Sistemas Dinamicos, Departamento de 
Matematica, Instituto Superior Tecnico, 1049-001 Lisboa, Portugal 



